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多 模式 集合 模拟 气候 变化 对 玉米 产量 的 影响 


摘 要 : 气候 模式 驱动 
差异 使 得 研究 
用 1981 一 2009 
作物 统计 模型 并 验证 了 
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作物 模型 是 气候 变化 影响 评估 的 主要 手段 .但 是 ， 单 一 气候 模式 输出 和 作物 模型 的 结构 
告 果 存在 不 确定 性 。 多 模式 集合 的 概率 预 估 可 以 有 效 减 少 研究 结果 的 不 确定 性 。 为 此 ， 本文 利 
年 东北 地 区 海伦 、 长 岭 、 本 溪 3 地 区 农业 气象 站 的 历史 气象 资料 和 玉米 作物 数据 ， 分 别 建 立 了 
APSIM 机 理 模 型 在 研究 区 域 的 适用 性 。 在 此 基础 上 ,与 CMIP5 在 RCP4.5 情景 下 的 8 


个 全 球 模式 结合 ， 尝 试 基于 多 模式 集合 评估 了 未 来 2010—2039 年 时 段 和 2040—2069 年 时 段 气候 变化 对 玉米 


量 的 可 能 影响 (相对 于 1976—2005 年 基准 时 段 )。 研 究 结果 表明 ，APSIM 模型 对 玉米 生长 发 育 和 产量 形成 有 


很 好 的 模拟 能 力 。 玉 米 生 育 期 的 模拟 误差 (RMSE) 为 3~4 d, 产量 的 RMSE 为 0.6~0.8 thm “。 建 立 的 产量 统计 
模型 表明 ， 玉米 出 音阶 段 (5 


旬 ) 的 温度 和 降水 的 增加 、 


光照 的 不 足 均 不 利于 产量 增加 。 与 1976—2005 年 基准 时 段 相 比 ,气候 因素 影响 1 


2010—2039 年 玉米 产量 减少 3.8%( 海 伦 )~7.4%( 本 溪 ), 减产 的 概率 为 64%( 长 岭 )~73%( 本 溪 ); 2040—2069 年 时 


月 中 旬 ) 的 温度 增加 对 产量 增加 有 积极 作用 ,而 开花 到 成 熟 阶 段 (7 月 中 旬 到 9 HJ 


7| Im 


段 减 产 6.499 (18 16)- 10.59 24 3X), Jà P^ US EE 29 74%( 海 伦 )~83%( 本 溪 )。 未 来 2010 一 2039 年 时 段 和 2040 一 2069 


HE 


时 段 基 于 统计 模型 得 


中 图 分 类 


基于 机 理 模型 模拟 的 产量 降低 分 别 为 6.6%( 海 伦 )~8.9%( 本 溪 ) 和 9.7%( 海 伦 )~13.7%( 本 溪 )， 均 高 于 相应 
到 的 0.9%( 海 伦 )~6.0%( 本 溪 ) 和 2.0%( 长 岭 )~7.3%( 本 溪 ) 减 产 结 
关键 词 : 气候 变化 ; 统计 模型 ; 机 理 模 型 ; 集合 模拟 ; 玉米 ; 生育 期 ; 产量 
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Abstract: Climate projections through process-based statistical crop models are important in studying the impacts of climate 


change on agricultural production. However, extensive assessments have generally relied on single climate with single crop 


models which have shown large discrepancies in predicted crop yields and estimations uncertainty hardly assessed. The proper 


understanding of uncertainties associated with such models is essential for effective use of projected results in devising adap- 


tation strategies. Assessing crop yield response to future climate conditions based on an ensemble of possible outcomes from 
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multiple climate projections and crop models could be more reliable than using a single model outcome. To estimate uncer- 
tainties associated with the study of the impacts of climate change on crop yield, we used 8 climate projections by GCMs un- 
der RCP4.5 in the CMIPS (which represented the uncertainties in the projected climate change) and a statistical process-based 
crop model (which represented the uncertainties in the different structures or different formulations of physiological processes 
of crop models). Historical data of crop and meteorological data during 1981—2009 from agro-meteorological stations of China 
Meteorological Administration in Hailun, Changling and Benxi in Northeast China were used to establish and evaluate statis- 
tical and process-based APSIM (Agricultural Production Systems sIMulator) models, respectively. Then the two crop models 
were linked with 8 climate projections to evaluate the impact of climate change on maize yield during 2010—2039 and 
2040-2069, using 1976-2005 as the baseline period. In total, 2 crop models under 8 climate projections for a period of 30 
years (a total of 480 simulations) were generated for both the baseline and two future climate periods. The results showed that 
APSIM model well simulated the growth and yield of maize. The root mean square error (RMSE) for the growth progress 
(flowering and maturity) simulation was 3-4 days and that for the yield simulation was 0.6—0.8 t-hm ?. The established statis- 
tical model suggested that temperature during emergence (mid May) had a positive effect on maize yield. However, the in- 
crease of temperature and rainfall, and lack of solar radiation during flowering and grain-filling periods (mid July to early 
September) had negative impact on increase of maize yield. Compared with 1976—2005, the resulting probability distributions 
indicated that due to climate change, maize yield in 2010—2039 decreased on average by 3.8% (Hailun) — 7.4% (Benxi), at a 
probability of 64% (Changling) — 73% (Benxi). During 2040-2069, maize yield increased by 6.4% (Hailun) — 10.596 (Benxi), 
at a probability of 74% (Hailun) — 83% (Benxi). The simulated yield decrease by the APSIM model was 6.6% (Hailun) — 8.996 
(Benxi) during 2010-2039 and 9.7% (Hailun) — 13.7% (Benxi) during 2040-2069. These were higher relative to those simu- 
lated by the statistical model, which were 0.9% (Hailun) — 6.0% (Benxi) during 2010-2039 and then 2.0% (Changling) — 7.396 
(Benxi) during 2040-2069. 

Keywords: Climate change; Statistical model; APSIM model; Ensemble simulation; Maize; Growth progress; Yield; 


作物 模型 是 评价 未 来 气候 变化 对 于 作物 产量 影 
响 的 主要 工具 。 其 中 ,机 理 模 型 因 具 有 较 强 的 机 理 


结果 。 如 ， 张建平 等 已 ] 应 用 WOFOST 模型 与 不 同 
气候 情景 数据 结合 模拟 东北 玉米 产量 时 得 到 ， 基 于 


性 ， 是 气候 变化 影响 研究 中 最 常用 的 工具 ， 如 
APSIMLII、DSSATP、WOFOSTBI、ORYZAD 等 。 然 
而 ， 其 复杂 的 过 程 需要 大 量 详细 的 作物 生长 、 管 理 、 
土壤 等 初始 数据 进行 模型 校 验 (1。 这 些 数 据 的 不 准 
确 性 , 增加 了 模拟 结果 的 不 确定 性 "中 ,相反 , 统计 模 
型 则 无 需 利 用 田间 生产 和 管理 数据 进行 模型 校 验 ， 
又 可 以 通过 有 限 的 数据 ， 潜 在 地 抓 住 相对 难以 理解 
的 某 些 过 程 (如 病虫害 的 动态 信息 ) 的 影响 , 有效 减 
少 了 机 理 模型 大 量 参 数 带 来 的 不 确定 性 *'" 相 对 于 
机 理 模型 在 气候 变化 影响 中 的 广泛 应 用 ,统计 模型 
的 作用 逐渐 受到 关注 。 例如, Lobel 等 中 建立 了 国家 
尺度 的 历史 气象 因子 和 玉米 (Zea mays)、 小 麦 
(Triticum aestivum), 7KTB(Oryza sativa). KE (Glycine 
max)4 种 作物 产量 的 统计 模型 ,量化 了 1980—2008 
年 气候 变化 对 全 球 粮食 作物 的 影响 。Tao 等 "通过 
统计 模型 分 析 了 1981—2009 年 气象 条 件 及 栽培 管 
理 措施 等 对 我 国 水 稳产 量 的 影响 及 贡献 ， 有 效 弥 补 
了 机 理 模型 量化 历史 气象 条 件 下 载 培 管理 措施 对 产 
量 贡献 的 不 足 。 

当前 气候 变化 影响 研究 中 ， 多 为 单一 气候 模式 
驱动 单一 作物 模型 得 出 的 结果 。 而 且 由 于 气候 模式 
和 作物 模型 结构 的 差异 得 出 不 同 甚至 相反 的 评估 


RegCM3 输出 的 多 情景 加 权 平 均值 简称 REA 情景 
未 来 2011 一 2050 时 段 玉 米 产量 减产 约 8%， 而 基于 
PRECIS 输出 的 A2 和 B2 情景 产量 分 别 增加 20% 和 
7%。 相 反 ， 袁 东 敏 等 5 同样 应 用 PRECIS 输出 的 
B2 情景 与 改进 后 的 东北 玉米 生长 模拟 模型 (结合 , 模 
拟 得 到 未 来 2011 一 2050 时 段 玉米 产量 以 减产 为 主 。 同 
FÉ, Guo 等 09 使 用 HadCM3 与 CERES 模型 结合 在 不 考 
E CO; 肥料 作用 下 ,研究 得 到 未 来 华北 小 麦 产量 增加 。 
而 Zhang 等 中 应 用 PRECIS 输出 驱动 APSIM 模型 模拟 
得 到 , 未 来 华北 地 区 小 麦 减产 。 这 些 问题 突出 了 在 气候 
变化 评估 中 考虑 由 气候 模式 和 作物 模型 的 结构 差异 引 
起 研究 结果 不 确定 性 的 重要 性 *”。 为 此 ， 本 文 以 东北 
地 区 黑龙 江 省 海伦 、 吉 林 省 长 肉 、 辽 宁 省 本 滩地 区 的 春 
玉米 为 例 ( 图 1), 借助 作物 机 理 模型 和 统计 模型 , 在 8 个 
GCMs 驱动 下 ,基于 多 模式 集合 模拟 了 未 来 2010 一 2039 
时 段 和 2040 一 2069 时 段 气候 变化 对 玉米 产量 的 可 能 影 
响 ， 旨 在 减少 单一 模型 对 模拟 结果 不 确定 性 的 影响 。 

1 材料 与 方法 

1.4. 研究 区 域 概况 及 数据 来 源 


东北 地 区 主要 包括 黑龙 江 、 吉 林 和 辽宁 三 省 。 
属 寒 温 带 湿润 、 半 湿润 气候 带 , 冬季 低温 干燥 , 夏季 温 
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暖 湿润 。 年 均 气 温 -3~10 C, 年 降水 量 400-1 000 mm, 
其 中 80% 集 中 发 生 在 5 一 9 月 。 土地 肥沃 ,有机质 含量 
达 12.4~34.2 gkg-!。 本 区 春玉 米 主要 是 一 年 1 熟 制 ， 
播种 面积 占 全 国 相 应 作物 总 播种 面积 的 30% 左 右 ， 
是 国内 重要 的 商品 粮 生产 基地 。 

本 研究 基于 能 够 代表 主要 玉米 产 区 的 生产 水 平 | 
和 县 有 较 完善 且 较 长 时 段 的 数据 资料 原则 ， 选 取 了 
黑龙 江 省 海伦 (47.4*N, 127.0°E)、 吉林 省 长 岭 (44.2°N,， 
124.0*E), U T $82 (41.3?N, 123.8"E)3 个 中 国 气 象 
局 农业 气象 试验 站 作为 研究 点 (图 1)。 研究 站 点 的 作 
物资 料 ， 来 源 于 国家 气象 信息 中 心 ， 包括 1981 一 45*N | 
2009 年 作物 生育 期 和 产量 ; 同期 历史 观测 气象 资料 
来 源 于 中 国 气象 局 科学 数据 共享 服务 网 包括 日 最 
高 温度 、 日 最 低温 度 、 日 照 时 数 、 降 雨量 。 其 中 , 日 
照 时 数 根据 Angstróm-Prescott 公式 转化 为 总 辐射 ?1。 
参照 玉米 的 常年 生长 发 育 期 进一步 求 取 生 育 期 内 jond 
逐 旬 的 平均 最 高 温度 、 平 均 最 低温 度 、 累 积 降 雨量 
和 累积 辐射 量 。 和 气候 情景 数据 来 源 于 国家 气候 中 心 ， 


55?N 4 


-z 


fu 
Changling 


-. 070140 280 420 560 
-———— — | 


选用 CMIP5 中 的 8 个 全 球 模式 BCC-CSM1-1( 中 国 )、 ds BUB EE 
CCSMA(ÉBI), CSIRO-MK3-6-0 GB ACRI NE). 图 1 研究 站 点 位 置 | 
EC-EARTH( 欧洲 ). GFDL-ESM2G( E js Fig.1 Locations of the selected study stations 


IPSL-CM5A-MR( 法 国 )、MRICGCM3( 日 本 ) 和 重播 值 的 方法 ， 降 尺度 到 研究 站 点 。 
NorESM1-M( 挪 威 )， 在 RCP4.5 情景 下 模拟 生成 未 来 12 研究 方法 

气候 时 段 (2010 一 2039 年 和 2040—2069 年 ) 和 基准 气 本 文 基于 收集 的 作物 、 气 象 等 数据 ,建立 作物 统 
候 时 段 (1976 一 2005 年 ) 逐 日 数据 ,包括 日 最 高 温度 、 计 模 型 并 验证 作物 机 理 模 型 ， 通 过 不 同 气候 模式 输出 
最 低温 度 、 降 水 。 模 式 分 辩 率 均 为 1x1。 辐 射 数据 ” ”驱动 不 同 作 物 模型 , 采用 集合 概率 方法 量化 气候 变化 
用 历史 观测 数据 代替 。 气 候 模式 数据 通过 反 距 离 权 ”对 玉米 产量 的 可 能 影响 。 具 体 技术 流程 见 图 2。 


历史 气象 和 作物 数据 
Historical meteorological and crop data 


8 组 气候 情景 数据 


8 climate projections data 


建立 作物 统计 模型 评价 作物 机 理 模型 
Establishment of crop statistical model Crop process-based evaluation model 


集合 产量 模拟 BC MS 
Ensemble yield CCSM4 
simulations CSIRO-MK3-6-0 
mE EC-EARTH 
概率 输出 GFDL-ESM2G 
Probabilistic output 
IPSL-CMSA-MR 
未 来 不 同时 段 玉 米 产 量 较 基 准时 段 变化 MRICGCM3 
Projected yield change relative to the baseline NorESMI-M 


2 气候 变化 对 玉米 产量 影响 研究 的 流程 图 
Fig.2 Flowchart for the study of climate change impacts on maize yield 
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1234 ”作物 机 理 模型 

APSIM 农业 生产 系统 模拟 模型 自 引 入 我 国 以 来 ， 
已 经 被 广泛 应 用 于 评价 农作物 生长 受气 候 波动 和 环境 
变化 的 影响 等 方面 中。 本 研究 基于 站 点 尺度 ， 以 单 
一 气候 区 的 生态 型 品种 作为 该 区 代表 性 品种 ， 应 用 
1981—1995 年 数据 ,采用 “ 试 错 法 ”， 对 模型 的 品种 
参数 进行 调试 , 使 模型 模拟 的 开花 期 、 成 熟 期 和 产 
量 能 够 最 大 程度 地 接近 观测 值 ， 应 用 1996 一 2009 年 
数据 对 校正 后 的 模型 进行 验证 。 选 择 以 下 统计 指标 
对 模型 的 模拟 能 力 进行 评价 ,包括 模拟 值 与 实测 值 
拟 合 方程 y-ax 的 决定 系数 R*(1 : 1D)， 反 映 了 模型 模 
拟 值 相对 于 实测 值 的 真正 偏差 ; 斜率 c， 反 映 了 模型 
整体 高 估 或 低估 实测 值 的 程度 ; 均 方 根 误差 RMSE, 
反映 了 模型 模拟 值 相 对 实测 值 的 绝对 误差 。 

基于 APSIM 模型 模拟 研究 站 点 代表 性 品种 气 
候 变化 下 的 产量 时 ， 播 种 期 采用 固定 播 期 ， 为 农业 
气象 试验 站 1981 一 2009 年 播种 期 的 平均 ， 收 获 期 为 
生理 成 熟 日 。 其 他 管理 措施 如 灌溉 和 施肥 均 设 置 为 
最 适 管理 ， 即 不 考虑 胁迫 影响 。 
1.2.2 ”作物 统计 模型 

应 用 1981 一 2009 年 玉米 产量 数据 建立 作物 统 
计 模 型 。 模 型 的 建立 包括 趋势 产量 (主要 受 农业 技术 


进步 等 的 影响 ) 和 气象 产量 (主要 受气 候 因 子 的 影响 )， 


以 便 和 机 理 模 型 输出 的 产量 同等 量 级 。 趋 势 产量 采 
用 多 项 式 拟 合 方法 得 到 。 


海伦 : Y, - 0.000 24? «0.007 067? +0.158 94t 2.34171 
R?—-0.83 (1) 
fll: y, -—0.000 837 0.03153? —0.183 01: +6.226 09 
R?-0.52 (2) 
Aki: Y, = —0.001 686 +0.072 446? —0.714 40t 4-6.091 59 
R=0.59 (3) 


式 中 : 7 是 趋势 产量 , :是 年 份 序列 ( 取 值 范围 为 1~29， 
分 别 代表 1981 年 至 2009 年 )。 由 于 趋势 产量 外 推 的 
限制 ， 本 研究 未 考虑 气候 变化 下 品种 、 栽 培 技 术 等 
的 更 新 对 产量 的 贡献 ， 等 于 最 近 2005 一 2009 年 趋势 
产量 平均 水 平 。 

从 作物 实际 产量 中 剔除 趋势 产量 ， 即 得 到 气象 
产量 。 气象 产量 基于 玉米 生育 期 内 64 个 旬 气 象 因子 
采用 逐步 回归 方法 拟 合 得 到 。 逐 步 方法 标准 使 用 F 
分 布 概率 ， 满 足 4=0.05 因子 引入 , 4=0.1 因子 剔除 。 
1.2.3 ”产量 分 析 

本 研究 使 用 8 个 气候 模式 数据 驱动 2 个 作物 模 
型 ， 在 基准 气候 1976—2005 年 时 段 、 未 来 2010— 
2039 年 和 2040 一 2069 年 时 段 ， 分 别 产 生 480 组 产量 
模拟 数据 。 考 虑 到 气候 变化 中 应 用 不 同 气候 模式 和 


不 同 作物 模型 对 研究 结果 的 可 能 影响 及 不 确定 性 ， 
研究 通过 多 模式 结合 ， 使 用 概率 密度 (PDFs，prob- 
ability density functions) 和 累积 概率 分 布 西数 (CDFs， 
cumulative distribution functions) 量 化 分 析 气 候 变化 
下 玉米 产量 相对 于 基准 情景 变化 。 
2 结果 与 分 析 
2.1 APSIM 模型 的 站 点 校 验 

3 显示 了 模型 模拟 的 玉米 开花 期 、 成 熟 期 及 
产量 与 观测 值 的 验证 结果 。 o 接近 1.0, 表明 模 型 模 
拟 值 与 观测 值 接近 。 生 育 期 和 产量 的 R >0.6， 表 明 
模型 可 以 解释 60% 以 上 的 观测 生育 期 和 产量 变化 。 

育 期 的 RMSE 为 3-4 d, 平均 误差 不 足 3%。 产 量 

的 RMSE 为 0.6~0.8 thm’, 平均 误差 小 于 11%。 整 
体 上 , 各 生育 期 间 的 RMSE 均 表 现 为 随 着 生育 进程 
的 增加 ,误差 增 大 。 这 主要 是 由 于 模型 模拟 的 生育 进 
程 取决 于 之 前 生育 进程 的 累积 。 而 生育 进程 的 模拟 误 
差 ， 也 会 影响 产量 的 模拟 误差 。 此 外 , 考虑 到 观测 过 
程 中 的 可 能 误差 ,如 取样 的 不 均匀 等 ,认为 APSIM 
模型 对 站 点 的 玉米 生长 具有 较 好 的 模拟 能 力 。 
22 ”作物 统计 模型 建立 

产量 统计 模型 中 的 各 参数 及 其 系数 如 表 1 所 
示 。 从 回归 模型 可 以 看 出 ， 玉 米 开花 灌浆 阶段 (7 月 
中 旬 到 9 月 上 旬 ) 的 温度 和 产量 具有 负 相 关 关 系 。 该 
阶段 温度 的 增长 加 速 了 生育 进程 , 减少 了 生物 量 积 
时， 从 而 影响 了 产量 形成 。 相 反 ， EKHE NEG 
月 中 旬 ) 的 温度 和 气象 产量 具有 正 相 关 关 系 。 该 阶段 
温度 的 适度 增长 利于 作物 出 苗 及 壮 苗 形 成 。 玉 米 开 
花 灌浆 过 程 中 大 量 降雨 或 光照 不 足 的 天 气 ， 容 易 影 
响 散 粉 和 授粉 ,诱发 籽粒 败 育 ， 最 终 导 致 产量 减少 。 
2.3 不同 作 物 模 型 模拟 的 玉米 产量 变化 

表 2 显示, 基于 不 同 作物 模型 模拟 的 玉米 产量 在 
未 来 气候 变化 下 均 降 低 。2040 一 2069 年 时 段 的 减产 
率 高 于 2010 一 2039 年 时 段 。 也 就 是 说 ， 随 着 气候 变 
化 的 加 剧 , 减产 率 增 大 。 基 于 机 理 模 型 模拟 的 未 来 两 
个 时 段 减产 率 均 高 于 统计 模型 的 模拟 结果 。2010 一 
2039 年 时 段 , 最 突出 的 为 长 岭 地 区 , 基于 APSIM 机 
理 模 型 的 产量 降低 约 8.4%， 而 统计 模型 的 模拟 结果 
{X 1.4%。2040 一 2069 年 时 段 ， 同样 以 长 岭 地 区 最 为 
明显 ， 基 于 APSIM 机 理 模型 的 产量 降低 约 12.6%, 而 
统计 模型 的 产量 降低 仅 为 2.0%。 从 研究 点 的 整体 平 
均 变 化 看 , 未 来 2010 一 2039 年 气候 变化 下 基于 机 理 模 
型 模拟 的 产量 降低 约 6.6%( 海 伦 )~8.9%( 本 溪 )，2040 一 
2069 年 时 段 产量 降低 为 9.7%( 海 伦 )~13.7%( 本 尝 ); 
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开花 期 Flowering 


220 


y= 0.99x 


210 R-073 R - 0.65 
RMSE-3.8 * RMSE-3.7 
200 y 0.99x 
R - 0.65 
RMSE=2.8 
190 1 1 1 
190 200 210 220 230 


成 熟 期 Maturity 


模拟 生育 天 数 Simulated days after sowing (d) 


5 " 1 i 
245 255 265 275 285 
观测 生育 天 数 Observed days after sowing (d) 


产量 Yield 


9r 
y= 1.08x 
R =0.61 
RMSE=0.75 

7 


R =0.63 
RMSE=0.80 


模拟 产量 Simulated yield (thm?) 


3 5 7 9 11 
观测 产量 Observed yield (thm?) 
图 3 海伦 (。)、 长 岭 (,) 和 本 溪 (e) 玉 米 开花 期 、 成 熟 期 天 


数 和 产量 的 观测 值 和 模拟 值 验证 结果 
Fig.3 Validation results between the observed and simulated 
days of flowering and maturity stages, and yield of maize at 
Hailun(®), Changling(4), Benxi(e) 


相反 ,统计 模型 的 模拟 产量 降低 程度 较 小 ， 相 应 时 


段 分 别 为 0.9%( 海 伦 )~6.0%( 本溪 ) 和 2.0%( 长 
岭 )~7.3%( 本 溪 )。 基于 不 同 模 型 模拟 产量 变化 的 标准 


偏差 及 离散 程度 没有 明显 差异 。 
24 气候 变化 对 玉米 产量 影响 的 概率 分 析 

基于 多 模式 集合 , 气候 变化 下 2010—2039 年 时 
段 和 2040—2069 年 时 段 玉 米 产量 相对 基准 时 段 变 
化 的 概率 密度 和 累积 概率 结果 (图 4) 可 看 出 , 产量 变 
化 概率 密度 符合 一 般 正 态 分 布 , 但 是 峰值 位 置 随 着 
气候 变化 的 加 剧 向 左 偏 移 增 大 ， 负 面积 增加 ， 表 明 
产量 降低 的 可 能 性 增 大 。 气候 变化 下 2010 一 2039 年 
时 段 玉 米 产量 降低 的 概率 约 为 64%( 长 肉 )~73%( 本 
31$), 2040— 2069 年 时 段 概 率 增 加 到 74%( 海 
伦 )~83%( 本 溪 )。 与 基准 时 段 玉米 产量 相 比 ，2010 一 
2039 年 时 段 玉米 产量 平均 降低 3.8%( 海 伦 )~7.4%( 本 
3€), 2040—2069 年 时 段 产量 继续 降低 ， 达 6.4%( 海 
伦 )~10.5%( 本 尝 )。 整 体 上 , 气候 变化 下 玉米 产量 变 
化 数据 离散 程度 较 高 , 平均 约 11%。 

3 讨论 与 结论 

不 同 结构 作物 模型 模拟 结果 表明 ， 基 于 机 理 模 
型 模拟 的 玉米 产量 变化 大 于 统计 模型 模拟 结果 。 这 
与 当前 的 相关 研究 结论 一 致 。 大 量 基 于 机 理 模 型 评 
估 的 产量 变化 以 -11%~20% 为 主 户 209， 而 统计 模型 
的 研究 结果 变化 不 足 10%。 如 马 玉 平等 中 "利用 积分 
回归 法 建立 的 我 国 不 同 省 区 玉米 产量 与 气象 因子 
的 相关 模型 ， 预测 未 来 全 国 玉米 产量 降低 在 7% 以 
Ao M Zhang SEU Zhou 等 2?] 应 用 CIMPS 的 多 
个 GCMs 数据 与 不 同 参数 的 作物 统计 模型 结合 ， 
评估 得 到 未 来 东北 地 区 玉米 产量 变化 不 足 5%。 此 
外 , Roudier 等 69 通过 综合 分 析 西 非 地 区 16 篇 气候 
变化 对 作物 产量 影响 的 文献 ， 也 得 到 了 相似 的 研 
究 结果 。 这 可 能 是 由 于 作物 机 理 模 型 相对 于 统计 
模型 而 言 更 能 捕获 气候 因子 之 间 的 相互 关系 (如 极 
端 气温 、 降 水 、 平 均 气 温 ) 导 致 的 复杂 的 非 线 性 产 
量 响应 。 

本 研究 对 未 来 玉米 产量 的 评估 ， 仅 考虑 了 单一 
的 气候 因素 的 影响 ,并 没有 考虑 随 着 气候 变化 作物 
栽培 措施 的 调整 以 及 CO, 浓度 增加 等 的 影响 。 已 有 
研究 表明 ， 随 着 气候 变化 ,采取 一 定 的 适应 性 措施 ， 
如 ， 品 种 的 更 替 5， 播 种 期 的 调整 [351， 均 有 利于 作 
物产 量 的 提高 。 此 外 , CO 浓度 升 高 利于 缓解 或 抵消 
温度 升 高 造成 的 负面 影响 ， 提 高 作物 产量 。 

多 模式 集合 是 降低 气候 变化 影响 评估 结果 不 确 
定性 的 有 效 方法 。 本 研究 从 气候 模式 结构 和 作物 模 
型 结构 差异 对 研究 结果 影响 不 确定 性 角度 出 发 ， 通 
过 多 模式 集合 评估 ， 在 一 定 程 度 上 反映 了 未 来 东北 
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R1 海伦 、 长 岭 和 本 溪 3 地 玉米 产量 统计 模型 参数 系数 


Table 1 Coefficients of parameters of maize yield statistical models at Hailun, Changling and Benxi 


气象 产量 Meteorological yield 


ns 趋势 产量 


ET Trend yield i Positive correlation fa 4H X Negative correlation 常数 项 
(thm ^) 参数 系数 参数 系数 Constant 
Parameter coefficient Parameter coefficient 
海伦 5 月 中 旬 最 低温 度 9 月 上 旬 最 低温 度 
Hailun 7.13 Average minimum temperature 0.009 Average minimum temperature in —0.209 2.370 
in middle May early September 
8 月 中 旬 累 积 降水 
Accumulated precipitation in middle -0.008 
August 
长 内 8 月 下 旬 最 低温 度 
Changling 7.86 Average minimum temperature in —0.344 5.378 
late August 
AXE 5 月 中 旬 最 高 温度 7 月 中 旬 最 高 温度 
Benxi 6.34 Average maximum temperature 0.238 Average maximum temperature in —0.290 0.801 
in middle May middle July 
8 中 旬 总 辐射 
Accumulated solar radiation in 0.014 


middle August 


表 2 基于 作物 统计 和 机 理 模 型 的 2010—2039 年 和 2040 一 2069 年 时 段 海伦 、 长 岭 、 本 溪 3 地 相对 于 基准 时 段 
(1976—2005 年 ) 的 5%、95% 和 平均 气候 变化 的 玉米 产量 变化 (Mean) 及 标准 偏差 (S.D.) 
Table 2 Maize yield changes (Mean) and standard deviation (S.D.) under 5%, 95%, mean climate changes during the periods of 
2010—2039 and 2040—2069, relative to 1976—2005 based on statistical and process-based crop model at Hailun, Changling, Benxi 


地 点 BER APSIM 机 理 模 型 Process-based APSIM model 统计 模型 Statistical model 
Site Period Mean (94) S.D.(96) P5 (96) P95 (%) ^ Mean (9*6) S.D. (94) P5 (94) P95 (96) 
海伦 2010 一 2039 —6.6 8.2 -24.9 3.0 -0.9 8.1 -13.1 8.9 
Hailun 2040— 2069 -9.7 10.6 -29.7 6.3 -3.1 7.6 -17.9 6.7 
Kus 2010— 2039 -84 10.2 -27.0 3.6 -14 15.2 -28.3 21.2 
Changling 2040— 2069 -12.6 10.0 -30.8 3.6 -2.0 7.2 -13.9 7.4 
X 2010— 2039 -8.9 12.1 -30.9 4.1 —6.0 11.9 -29.7 13.4 
Benxi 2040— 2069 -13.7 10.5 -32.0 3.8 -7.3 11.1 -27.6 8.2 


06 f 海伦 Hailun 2010-2039 1.0 
海伦 Hailun 2040-2069 
一 一 长 岭 Changling 2010-2039 
---- 长 岭 Changling 2040-2069 
一 一 本 溪 Benxi 2010-2039 
-~--- 本 溪 Benxi 2040-2069 


0.8 


累积 概率 CD 
o 
ON 


e 
P 


0.2 


产量 变化 Yield change (%) 


图 4 海伦 、 长 岭 、 本 溪 2010 一 2039( 实 线 ) 和 2040 一 2069( 虚 线 ) 时 段 玉 米 产量 相对 基准 时 段 (1976 一 200s 年 ) 变 化 的 概率 
密度 及 累积 概率 分 布 


Fig.4 Probability density (PD) and cumulative density (CD) for simulated maize yield changes during 2010-2039 (solid line) and 
2040—2069 (dash line), relative to baseline (1976—2005) at Hailun, Changling, Benxi 
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地 区 玉米 产量 变化 的 可 能 范围 ， 得 到 未 来 2010 一 
2039 年 时 段 和 2040 一 2069 年 时 段 玉米 产量 分 别 降 
低 3.8%~7.4% 和 6.4%~10.5%。 这 与 Tao 等 9 基于 
多 套 气候 情景 数据 与 作物 机 理 模 型 结合 的 研究 结果 
相似 。 值 得 注意 的 是 ，Zhou 等 中 通过 综合 考虑 气候 
模式 和 作物 统计 模型 内 部 参数 引起 的 不 确定 性 ， 评 
估 得 到 未 来 东北 地 区 玉米 产量 有 增加 趋势 。 这 些 研 
究 差异 体现 了 作物 模型 选择 的 敏感 性 ， 另 一 方面 也 
突出 了 考虑 模型 内 部 参数 的 重要 性 。 因 此 ， 结 合 前 
人 研究 成 果 ， 综 合 考虑 气候 模式 结构 、 作 物 模型 结 
构 及 模型 内 部 参数 不 确定 性 对 研究 结果 的 影响 ， 是 
下 一 步 研究 的 重点 。 


致谢 ， 户 心 感谢 国家 气候 中 心 徐 颖 研究 员 为 本 研究 
提供 的 气候 情景 数据 支持 。 
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